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ABSTRACT 

Extensive photoionization model grids are computed for single star H II regions using stellar atmosphere models 
from the WM-Basic code. Mid-IR emission line intensities are predicted and diagnostic diagrams of [Ne III/ II] 
and [S IV/III] excitation ratio are build, taking into account the metallicities of both the star and the H II region. 
The diagrams are used in conjunction with galactic H II region observations obtained with the ISO Observatory 
to determine the effective temperature r e ff of the exciting O stars and the mean ionization parameter U. T e {{ and 
U are found to increase and decrease, respectively, with the metallicity of the H II region represented by the 
Ne/Ne ratio. No evidence is found for gradients of T e ff or U with galactocentric distance R ga i. The observed 
excitation sequence with R ga i is mainly due to the effect of the metallicity gradient on the spectral ionizing 
shape, upon which the effect of an increase in r e g with Z is superimposed. We show that not taking properly 
into account the effect of metallicity on the ionizing shape of the stellar atmosphere would lead to an apparent 
decrease of T e ff with Z and an increase of r e ff with R ga i. 

Subject headings: Galaxy: stellar content — infrared: ISM — (ISM:) HII regions — stars: atmospheres — 
stars: fundamental parameters — (stars:) supergiants 



1. INTRODUCTION 

The determination of the stellar distribution (especially of 
the hottest stars) and physical characteristics of galactic H II 
regions are of primary importance to evaluate star formation 
theories and for our understanding of the chemical evolution 
o f galaxies. 

IShields & Tinslevl (119761) used the equivalent width of the 
Hg emission from H II regions in spiral galaxies to determine 
the existence of a radial gradient in the effective temperature 
(hereafter T e ff) of the hottest stars , assoc iated with a decrease 
in metal abundance. Campbell (1988) determined T e ff and 
the ionization parameter U for various H II galaxies, and 
concluded that the T e ff of the hottest star decreases with 
increasing oxygen ab undance. 

On the other hand, Evans & Dopita ( 1985) have computed 
extens ive photoionization models, using Hummer & Mihalas 
( 1970) atmosphere models, and have determined from optical 
observations of H II regions that the ionization temperature 
of the exciting stars is approximately constant (41 kK, 
independently of the metallicity Z and U). They found, 
howev er, an anticorrelation between U and Z. Fierro et al. 
(1986) found a near constant T & g of 35 kK between 1 and 
5 kpc from the cent er of NGC 2403. 

More recently, Mart in-Hernandez et alJ ( l2002bl) used 
Infrared Space Observatory (ISO) spectral observations of 
galactic H II regions to show that the gas excitation increases 
with the galactocentric distance R ga i. They concluded that 
the stellar spectral energy distributions (hereafter SEDs) are 
softer at higher metallicities, that is towards the galactic center 
and that the SEP cha nges can explain the observed gradient. 
iGiveon et alJ (|2002b) similarly used ISO observations, but 
suggest that the increase in excitation corres pond instead to 
a decr ease in stellar effective temperature. iMorisset et al.l 
d2003l hereafter MSBM03) show that excitation gradients are 
partly due to changes in the ionizing spectral shape of O stars 
with metallicity. They concluded that the excitation scatter is 
probably mainly due to randomization of both the stellar r e ff 



and the nebular mean ionization parameter U . 

No attempt was ma de by IGiveon et alJ 

d2002bl) . iMartm-Hernandez et all d2002bl) . nor in MSBM03 
to determin e 7ht and U f or individual H II regions. 

iDors & Copettilll200 3) used optical observations of galactic 
and magellanic H II regions to determine T e ff from optical 
diagnostic line ratios. They also found an increase of T e ff with 
the galactocentric distance. 

The aim of the present work firstly is to build diagnostics 
diagrams for the determination of T e s and U, based on mid-IR 
emission lines. The diagrams are derived from a extensive 
grid of photoionization mod els that populate the T e s-U-Z 
space and use the WM-Basic (Pauldrachet al. 2001) code to 
compute the ionizing atmosphere models. In a second step, 
T e ff and U are determined for the ISO H II regions using the 
new diagnostic diagrams. 

Sect. |2 describes the ISO observations of H II regions, and 
Sect. 13 the grid of photoionization models. The location 
of ISO observations in the model grids, and the process to 
determine T e g, and the mean ionization parameter U for every 
object are presented in Sect.|4l using two different methods. 
Sect. |3 describes the resulting gradients of T e s and U. The 
effect of the stellar metallicity in the determination of J e ff is 
discussed in Sect. [6] in particular the influence of the changes 
in the stellar SEDs with metallicity. The conclusions are 
presented in Sect.Q 

2. ISO OBSERVATIONS OF H II REGIONS 

Mid-IR fine-structure line intensities obtai ned from 
observ a tions of H II regions w i th ISO -SWS ( IGiveon et alJ 
2002b; Martin-Hernandez et al. 2002a) are used here to 
determine the various stellar and nebular parameters. 

Line ratios of [ArIII]8.98^m/ [ArII]6.98/im, [SlV]10.5^m 

1 T he mean ionization parameter U is defined following Evans & Dopita 
i 1985) as the value of U evaluated at a distance from the ionizing star 
r = r emply + AR/2, where r empty is the size of the empty cavity and AR is 
the thickness of the uniform density H II shell. 
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/ [S III] 18.7/Ltm, and [Ne III] 15.5/im/ [Ne II] 12.8/im (hereafter 
[Ar III/II], [S IV/III] and [Ne III/II] resp.) ai-e then used to build 
excitation diagrams (see MSBM03 for more details). The 
line int ensities were corrected for reddening bv lGiveon et alJ 
(2002a). Once the sources for which at least one of the line 
intensities used in this work is not defined (or have only an 
upper limit) are removed, 42 usable sources remain. 

3. GRID OF PHOTOIONIZATION MODELS 

A grid of photoion ization models using the NEBU code 
( Morisset et al. 2002) has been calculated using ionizing 
spectral distributions from supergiant atmosp here models that 
were computed with WM-Basic V. 2. II 2 tPauldrach et alJ 
1200 11) . The current grid of photoionisation models is similar 
to the one used in MSBM03 (see MSBM03 for details). It 
has been extended further to encompass the full range of 
values expected within the 3D parameters space T^f-U-Z. T e s 
is ranging from 30 to 50 kK, by 5 kK steps, log(t7) being 
-2.6, -1.5, -0.8, -0.1 and 0.5, and the metallicity Z of both 
the ionizing star and the nebular gas being 0.5, 0.75, 1.0, 1.5 
and 2.0 times the solar value (as defined in WM-Basic). In 
total, 125 models have been computed from which mid-IR 
line intensities were derived. 

4. DETERMINATION OF r cff AND U 

Three mid-IR line ratios could in principle be used as 
excitation diagnostics, namely [Ar III/ II], [Ne III/ II] and 
[S IV/III]). We note however that [Ar III/II] is overestimated in 
photoionization models, as pointed out in MSBM03. It is not 
clear whether the latter is due to an improper determination of 
the ionizing flux near 24 eV or simply to incomplete atomic 
physics used within photoionization codes (e.g. missing 
accurate dielectronic recombination rates). For these reasons, 
despite its low dependence on U, the [Ar III/II] ratio proves to 
be useless for determining T e g. On the other hand, while the 
alternative excitation ratios [Ne III/II] and [S IV/III] are both 
sensitive to T e ff and U, it turns out that the effect of either 
parameter on both ratios is somewhat different, and a way to 
determine T e s and U using these excitation diagnostics can 
be extracted. Even though the direct use of [Ne III/II] and 
[S IV/III] or other combinations of these line ratios would 
provide equivalent constrain ts; we pre fer to adopt %-Ne> 
defined following Vilchez & Pagel ( 1988) as ?y s _ Ne = [S IV/III] 
/ [Ne III/II], in combination with [Ne III/II]. The determination 
of T e ff and U turns out to be clearer and easier to read off when 
using ?7s-Ne- 

The underlying hypothesis/assumptions made here are the 
following: 1) the ISO H II regions are excited by a single 
star (i.e. the presence of other less luminous stars doesn't 
affect the results), 2) the H II regions are ionized by stars 
of comparable surface gravity log(g) at a given r e g- (see 
MSBM03 for the effects of log(g)), 3) the Ne abundance 
determined from an H II region is a reliable estimator of 
the ionizing star metallicity, 4) the presence of dust in the 
H II regions doesn't affect the T e g and U gradients. Dust 
would decrease the global amount of ionizing photons, but 
also increase the hardness of the ionizing SED, increasing 
the excitation of the gas depicted by the IR excitation 
diagnostics (see MSBM03) and therefore the T e ff we infer can 
be overestimated, 5) Using a nebular geometry consisting of 
a simple shell does not affect the diagnostic diagrams, 6) the 
WM-Basic atmosphere models describe well the ionizing flux 

2 freeware available at http://www.usm.uni-muenchen.de/people/adi/ 
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FIG. 1. — The rjs-Ne versus [Ne III/II] excitation diagram for the H II 
regions observed by ISO. The grid of photoionization model with stellar and 
nebular metallicity given by Z/Zq = 1.0 (S-Method)_is drawn. Solid (dotted) 
lines connect iso-j7 (iso-r c ff) models. T e g and log(£7) take values of 30, 35, 
40, 45, 50 kK and -2.6, -1.5, -0.8, -0.1, 0.5, respectively. 



between 35 and 41 eV (or at least the relative changes that 
occurs when the parameters T e a or Z are varied). 

4.1. S-Method : Using only solar metallicity atmosphere 

models 

In a first step, we use only the results of the photoionization 
models obtained with the solar abundances atmosphere 
models. 

The Fig.^shows the values taken by ?/s-Ne and [Ne III/II], 
when r e ff and U are varied in locked steps. For each H II 
region, 2D-interpolations within this grid are performed and 
used to determine r e ff and U. All the observed values lie 
inside the grids, no extrapolation is needed. The T e s and 
U obtained with solar metallicity atmosphere models are 
presented in Sec.|5] 

From Fig. [2 we can determine the effects of uncertainties in 
line intensities on T e s andf/ : a factor of two in the excitation 
diagnostics leads to a shift in T e g by 1 kK and in U by 0.5 dex. 

4.2. Z-Method : T e ff and U obtained using Z-dependent 

atmosphere models 

The chemical composition of a star strongly affects 
its radiation, especially in the EUV, where the ionizing 
photons are emitted. For the same r e ff, changing the 
stellar luminosity by a factor of 4 can affect the excitation 
diagnostic line ratios by up to 2 orders of magnitude 
(MSBM03). The determination of r eff and U need then to 
be performed using grids of photoionization models with 
various stellar metallicities, as described in Sect. [3] Making 
the assumption that the metallicity of the ionized region 
reflects the metallicity of the ionizing star, we can use 
adapted diagnostic grids, corresponding to the H II regions 
metallicities, to determine r e ff and U. The metallicity of an 
H II region is hereafter given by the abundance ratio Ne/Ne^ = 
[Ne/H]/[Ne/H] Q , where [Ne/H] is obtained from lGiveon et alJ 
(2002a). The solar abundance used, [Ne/H] Q = 1.4 x 10" 4 , is 
defined a s the value of the abu ndance gradient [Ne/H](R ga i) 
found by iGiveon et alJ ( l2002al) . evaluated at 8.5 kpc. The 
Ne abundances as a measure of Z/Z Q , are preferred to a 
combination of Ne, Ar and S abundances, since, for these 
two last elements, the abundances are not reliable when the 
excitation is extreme (MSBM03). For each H II region, 
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we extract from the T e ff-U -Z photoionization model grid the 
Teff-U plane, whose metallicity lies closest to the Ne/Ne© of 
the H II region. The Fig.[2]shows, for the 5 metallicities used 
in the photoionization model grid, the values taken by r;s-Ne 
and [Ne III/ II]. The effect of increasing the metallicity of the 
atmosphere models (and consequently of the nebular gas) is 
clearly to increase the value of T e g for a given [Ne III/ II] 
ratio, as already pointed out in MSBM03. The observed 
values for the H II regions are also plotted in the diagramm 
corresponding to their metallicities. All the observed values 
lie inside the grids. 2D-Interpolations are performed to 
determine T e s and U for each H II region. 

5. RESULTS 

The Table ^ presents the characteristics of the 42 H II 
regions used in this work: R ga ], Ne/Ne0 (with a sy mbo l 
corresponding to the grid used within the Z-Method, Sec. l4.2K 
[Ne III/ II], [S IV/ III], and the T et£ and obtained using the 
Z-Method. The T eS range from 34 to 50 kK, and log(I7) 
from -1.5 to 0.5, with mean values of 40.5 kK and -0.60, 
respectively. Such range for T e a and U is in agreement with 
the results found by Evans & Dopita (119851) . For the 3 sources 
for which two independent observations are available, the 
results obtained lead to a coherent determination of T e s, while 
the values of U can differ by a factor up to 5. 

The set of inferred values for r e ff and log(f?) versus the 
galactocentric radius R„ a i and the abundan ce ra tio Ne/Ne© are 
shown in Fig. |3] for the S-Method (Sec. PTTV and in Fig. g] 
for the Z-Method (Sec. 14. 2> . Linear fit to the data are also 
presented in all the figures. 

The results obtained for T e ff with S- a nd Z-Methods are 
very different (upper panels of Figs. I3I4> . While S-Method 
leads to an increase (decrease) of T e g with R ga i (Ne/Ne Q ), 
the use of the Z-Method leads to rather different results: 
no clear correlation is found to exist between r e ff and R„ a i 
while a correlation is present between T e ff and Ne/Ne Q . 
The distribution of 7^ versus Ne/Ne Q obtained with the 
Z-Method can also be described as follows: T e g increasing 
strongly with Ne/Ne Q for Ne/Ne Q < 1 .2, and a quasi constant 
value (~ 43 kK) for higher Ne/Ne Q , coupled with a higher 
dispersion. 

The dispersion of Ne/Ne Q with position in the Galaxy 
is relatively high (see the symbols dispersion along R„ a i 
in the left panels of Fig. 0}. The absence of correlation 
between T e s and R ga i obtained with Z-Method might be the 
result of this high dispersion. Note that for H II regions 
with R„ a ]<R„ a ]^, the d etermination of R ga i is degenerate 
(Martin-Hernandez et al. 2002a) and the errors are also 
important. 

Virtually no changes are obser ved f or log(£/) from S- to 
Z-Method (lower panels of Figs. 1314-i . This gives insights 
of the robustness of our results concerning U whatever the 
method used. This can be understood as follows: the main 
effect of changing Z on the diagnostic diagrams (Fig. [2} is to 
shift horizontally ([Ne III/ II] axis) the grids of models, while 
the determination of U is mainly dependent on the vertical 
position along the f?s-Ne axis, which is not affected by the 
stellar metallicity. No clear correlation is found between U 
and R ga i but an inverse correlation between U and Ne/Ne© is 
present (see Fig.|4] lower panels). 

Fig. |5] shows the distribution of T e ff versus the [Ne III/ II] 
excitation ratio. The softening of the stellar emission when 
the metallicity increases (even if T e ff doesn't change) is 
enough to cover the whole observed range of [Ne III/ II] (over 
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Name" 




Ne/NeQ 1, 


[Ne III/ II] 1 ' 


[S IV/III] rf 


T eS eJ 


Of 


GCRINGSW 


0.00 


2.221= 


0.048 


0.021 


39.5 


0.302 


ARCHFILNW 


0.00 


2.340= 


0.037 


0.016 


38.8 


0.280 


IRAS 17455-2800 


0.50 


1.645° 


0.405 


0.156 


42.3 


0.213 


SGR D HII 


1.50 


1.410° 


0.398 


0.161 


42.2 


0.241 


WBH9815567-5236 


4.30 


0.773 w 


0.355 


0.343 


35.5 


2.521 


RAFGL 2094 


4.40 


1.647° 


0.056 


0.023 


36.1 


0.297 


IRAS 18317-0757 


4.50 


1.895= 


0.060 


0.027 


40.1 


0.314 


WBH9818317-0757 


4.50 


2.139= 


0.061 


0.020 


40.2 


0.163 


GAL 033.91+00.11 


4.50 


0.970+ 


0.482 


0.187 


38.2 


0.282 


IRAS 15502-5302 


4.60 


0.932+ 


0.216 


0.126 


36.3 


0.630 


IRAS 18434-0242 


4.60 


1.755= 


0.251 


0.165 


45.2 


0.508 


IRAS 18502+0051 


4.80 


1.987= 


0.204 


0.127 


44.5 


0.437 


GRS 328.30+00.43 


4.80 


0.741 w 


0.239 


0.177 


34.9 


1.148 


IRAS 17221-3619 


5.20 


1.706° 


0.036 


0.036 


34.9 


2.007 


WBH9816172-5028 


5.60 


1.707° 


0.327 


0.132 


41.5 


0.244 


GAL 337.9-00.5 


5.80 


2.076= 


0.578 


0.147 


48.2 


0.060 


G327.3-0.5 


6.30 


1.912= 


0.405 


0.078 


47.3 


0.038 


WBH9817059-4132 


6.30 


2.386= 


0.237 


0.039 


45.5 


0.034 


GAL 045.45+00.06 


6.30 


1.450° 


0.905 


0.381 


45.3 


0.268 


IRAS 15384-5348 


6.40 


1.384° 


0.296 


0.168 


41.0 


0.500 


GRS 326.44+00.91 


6.50 


1.284° 


0.308 


0.152 


41.2 


0.372 


WBH98 15408-5356 


6.60 


1.768= 


0.695 


0.189 


48.7 


0.066 


M17IRAMPOS8 


6.80 


1.363° 


2.839 


0.646 


49.7 


0.074 


W51 IRS2 


7.30 


1.106+ 


1.743 


0.807 


41.1 


0.596 


NGC6357I 


7.70 


1.169+ 


1.107 


0.407 


40.1 


0.270 


S106 IRS4 


8.40 


0.742 w 


0.115 


0.099 


34.0 


1.326 


NGC3603 


8.90 


1.524° 


2.513 


0.873 


48.3 


0.240 


IRAS 12063-6259 


9.30 


0.775 w 


1.930 


1.071 


38.4 


1.242 


GAL 289.88-00.79 


9.30 


0.765 w 


0.267 


0.159 


35.0 


0.716 


IRAS 10589-6034 


9.50 


0.853 w 


0.265 


0.165 


35.0 


0.799 


RAFGL 4127 


9.60 


0.827 w 


4.116 


1.654 


39.7 


0.759 


BE83IR 070.29+0 


9.60 


0.547 x 


1.954 


0.721 


34.4 


0.377 


BE83IR 070.29+0 


9.60 


0.414 x 


2.115 


0.566 


35.1 


0.086 


IRAS 11143-6113 


9.70 


1.037+ 


4.571 


1.226 


43.7 


0.342 


IRAS 19598+3324 


9.80 


0.432 x 


2.192 


0.784 


34.5 


0.329 


IRAS 12073-6233 


10.10 


0.733 w 


5.813 


2.368 


40.2 


1.093 


GAL 298.23-00.33 


10.10 


0.627 w 


6.110 


2.091 


40.4 


0.803 


IRAS 02219+6152 


11.00 


0.986+ 


2.402 


0.755 


42.1 


0.317 


IRAS 02219+6152 


11.00 


1.016+ 


2.252 


1.182 


41.6 


0.887 


W3 IRS2 


11.30 


0.961+ 


2.339 


0.665 


42.1 


0.250 


W3 IRS2 


11.30 


0.943+ 


2.321 


1.139 


41.7 


0.820 


SI 56 A 


11.50 


0.820 w 


0.161 


0.094 


34.4 


0.607 



Table 1 

Sources parameters and corresponding T eS and U 

D ETERMINED IN THIS PAPER. IRAS N AMED SOURCES FROM 

IMARTIN-HERNAnDEZ ET Al1 ( |2002aD . OTHERWISE FROM 
IGIVEON ET AL.H2002bI) . °: R„ a i IN KPC. c : UPPER SYMBOLS 
INDICAT WHICH GRID IS USED TO DETERMINE T e ff AND U: X, W, 
+ , O, = CORRESPOND TO THE NE/NE q = 0.5, 0.75, 1.0, 1.5 
AND 2.0 GRID RESP. (SEE FlG.|2j * : LINE INTENSITIES 
CORRECTED FOR REDDENING, SEE TEXT. ' : T cff IN KK. 1 : 

Values obtained using the Z-Method (SEC. l4.2t . 



2.5 dex). This Fig.|5]illustrates one more time the illusion of 
determining T e s from only one excitation diagnostic. On the 
other hand, a global increase of T e g with [Ne III/ II] is also 
present (the hottest stars are associated with the most excited 
H II regions). 

6. DISCUSSION 

The results presented in Sect.|5]are sensitive to the changes 
of the stellar SED with metallicity, for a given atmosphere 
model code. The gradients in T e g as a function of Ne/Ne Q 
and R ga i derived from a single solar metallicity set of stars 



4 



Morisset, C. 




FIG. 2. — Same as Fig.^but the stellar and nebular metallicities beeing 0.5, 0.75, 1.0, 1.5, 2.0 times solar, from upper left to lower right. Solid (dotted) lines 
connect iso-U (iso-T e ff) models. T c [{ and log(J/) take values of 30, 35, 40, 45, 50 kK and -2.6, -1.5, -0.8, -0.1, 0.5, respectively. All the observations of the H II 
regions are distributed in the diagramms, according to their nearest values of Ne/NeQ (Z-Mefhod). The H II regions are symbolized using X, W, +, 0, and = for 
Ne/NeQ = 0.5, 0.75, 1.0, 1.5, 2.0, respecively. See Tab.fTlfor the correspondance between sources and the diagramm used. 



(S-Method) are very different from the gradients obtained 
with coherent metallicity for the stellar atmosphere model 
(Z-Method). 

The results obt aine d with the S-M ethod agree with those of 
Campbell dl988l) and iDors & Copettil d2003l) . who similarly 
didn't take into account the Z-dependence of the stellar 
emission. This confirms that the trends seen in the 
upper panels of Fig. 0] when the Z-dependence is properly 
considered, are mainly due to the changes in the stellar SEDs 
when the metallicity decreases (an effect which has to be 
present, but whose magnitude might depend on the family 
of atmosphere models used). Dors & Copetti (2003) check 
the effect of metallicity on the excitation of the ionized gas, 



but didn't found a strong effect; the metallicity range they 
test is only a factor of two for WM-Basic models, they used a 
low value for the ionization parameter, the effect of Z being 
reduced in this case, and they used dwarf WM-Basic models, 
for which the effect of Z is lower than for supergiant s. The 
T e ff va lues obtained by Campbell (1988) and lDors & Copettil 
(2003) are derived using a S-Method, and the apparent 
decrease of T e ff with Z and increase of T e s with R ea i that they 
respectively found is likely to be the result of not considering 
Z-dependent stellar models (the excitation of the gas is not 
a valid T e ff indicator, as shown by Fig.[5}- Note also that the 
maximum T e s obtained with S-Method is 45 kK, while the T e ff 
values from Z-Method get as high as 50 kK. 



Effective temperature of galactic O Stars 



5 




FIG. 3. — r c ff (upper panels) and U (lower panels) versus R 6a i (left) and Ne/Ne0 (right) for the sample of H II regions. The grid of photoionization models 
results computed with Z = Zq stars show in Fig.Qis used for every source (S-Method, Sec. 14. II . Linear fits to the data are presented. 




FIG. 4. — Same as in Fig. [3] but each H II region is symbolized accordi ng to its abundance Ne/NeQ as in Fig. [2] The determination of 7" c ff and U follows from 
the use of the appropriate diagnostic diagram from Fig.|5](Z-]VIethod, Sec. 14.21 . Linear fits to the data are presented. 
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2.5x10 

0.001 0.010 0.100 1.000 10.000 100.000 
[Nolll]/[Ncll] (<1 1 .1 eV) 

FIG. 5. — r c ff versus [Ne III/ II] for the ISO H II regions. Symbols as in 
FigEI 



The results shown in this paper concerning T e g and U are 
strongly dependent on the kind of atmosphere model used to 
compute the photoionization grid of models. In MSBM03, we 
discuss in detail the major differences between, f or instance, 
WM-Basic and CMFGEN JHillier & Miller! 1998t) . in sofar as 
the determination of r en - is concerned. Using the CMFGEN 
instead of WM-Basic atmosphere models would certainly lead 
to a globally lower 77 e ff (MSBM03), and perhaps different 
gradients. 

The increase of 77 e ff with Z found h ere is in contra diction 
with the theoretical predictions of e.g.[Schaller et al. ( 1992), 
and if confirmed it might have profond implications for the 
study of the upper mass end of the IMF. 

However, more extensive studies will be needed to check 
whether the use of different atmospheres codes will confirm 



the gradients found here or generate genuine gradients. Our 
results, nevertheless, show the importance of taking properly 
into account the variation in the stellar SEDs with metallicity 
in any attempt to determine a reliable T e s from H II regions. 



7. SUMMARY AND CONCLUSION 

Based on WM-Basic atmosphere models we have computed 
a large set of photoionisation models. From these models 
we have built excitation diagnostic diagrams based on 
[Ne III/II] and [S IV/III] (mid-IR lines) excitation ratios. ISO 
observations of galactic H II regions are superimposed to 
these diagrams. According to their metallicity, 77 e ff and U are 
determined for every H II region. 

A correlation between T e ff and Ne/Ne Q , and 
an anti-correlation between U and Ne/Ne Q , have been found, 
without evidence of any correlation between both 77 e ff and U 
versus R ga i. The determination of T e s is strongly dependent 
on the changes in stellar SEDs due to the radial metallicity 
gradient within the Galaxy, while the results found concerning 
the behaviour of U are globally insensitive to this effect. The 
gaseous excitation sequence is therefore mainly driven by the 
effects of metallicity on the stellar SEDs. A global increase 
of r e ff with metallicity appears nevertheless to be present. 
More investigation using different atmosphere codes will be 
needed to confirm that our conclusions are not unduly biased 
toward the use of WM-Basic models. Comparison with T e ff 
determined from direct observations of ionizing stars can also 
help to evaluate the robustness of the method presented in this 
work. 
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